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ABSTRACT 

Mass-to-light versus colour relations (MLCRs), derived from stellar population 
synthesis models, are widely used to estimate galaxy stellar masses {Mf) yet a detailed 
investigation of their inherent biases and limitations is still lacking. We quantify sev¬ 
eral potential sources of uncertainty, using optical and near-infrared (NIR) photometry 
for a representative sample of nearby galaxies from the Virgo cluster. Our method for 
combining multi-band photometry with MLCRs yields robust stellar masses, while 
errors in decrease as more bands are simultaneously considered. The prior as¬ 
sumptions in one’s stellar population modelling dominate the error budget, creating 
a colour-dependent bias of up to 0.6 dex if NIR fluxes are used (0.3 dex otherwise). 
This matches the systematic errors associated with the method of spectral energy dis¬ 
tribution (SED) fitting, indicating that MLCRs do not suffer from much additional 
bias. Moreover, MLCRs and SED fitting yield similar degrees of random error (~0.1- 
0.14 dex) when applied to mock galaxies and, on average, equivalent masses for real 
galaxies with M* ~ 10*“^^ M©. The use of integrated photometry introduces addi¬ 
tional uncertainty in M* measurements, at the level of 0.05-0.07 dex. We argue that 
using MLCRs, instead of time-consuming SED fits, is justified in cases with complex 
model parameter spaces (involving, for instance, multi-parameter star formation histo¬ 
ries) and/or for large datasets. Spatially-resolved methods for measuring M* should be 
applied for small sample sizes and/or when accuracies less than 0.1 dex are required. 

An Appendix provides our MLCR transformations for ten colour permutations of the 
grizH filter set. 

Subject headings: galaxies: clusters: individual: Virgo; galaxies: general; galaxies: 
formation; galaxies: fundamental parameters; galaxies: stellar content 



- 2 - 


1. Introduction 


Stellar mass (M*) is one of the most fundamental parameters for tracing the formation of 
galaxies, as embodied in its many correlations with other galaxy properties, such as halo mass. 


size, star formation rate, gas-phase and stellar metallicity, at all redshifts (Tremonti et al. 2004 


Gallazzi et al.||2005[|Behroozi et al.||2010[|Mannucci et alT]|2010| ). Knowledge of M* also informs 


us on a variety of phenomena pertaining to galaxy formation efficiencies, dynamics, and feedback 
histories. Given its astrophysical relevance, there is natural interest in measuring M* as accurately 
as possible, whilst also achieving a complete account of its associated errors. 


For local galaxies, whose distances are fairly well-known, determining M* essentially con¬ 
sists of measuring the mean mass-to-light ratio (M*/L) of their stellar populations. The standard 
approach of measuring A/*/L is to fit stellar population synthesis (SPS) models to the SED, line 
strength indices, and/or full spectrum of a galaxy. Alternatively, one may access M^/L through 
maximal-disk fits (e.g. van Albada & Sancisi 1986|, relations between M^/L and broadband colour 


(e.g. [Bell et ali]|2003| ), or a combination of the two ( |Bell & de Jong||2001| ). These methods allow 
one to access M*/L whenever multi-band photometry or high-quality spectroscopy is unavailable. 

In addition to minimal input, M*/L-colour relations (MLCRs) make determining mass-to- 
light ratios quite trivial, through the use of (quasi-)linear relations. The linearity of MLCRs stems 
from the fact that changes in age, metallicity, and reddening scatter model stellar populations along 
the relations rather than away from them; equivalently, the age-metallicity-reddening degeneracy 
conspires to keep the scatter in MLCRs low (Bell & de Jong|[2001|). The simplicity of MLCRs 


has precipitated their widespread use, as well as several calibrations of these relations (Bell & de 


JongpOOl] |Bell et all]|2003t [Zibetti et'aL]|2009[ [Gallazzi & Bell|[2009| [Taylor et'aL]|2011[ [Into & 
Portinari|2013[|McGaugh & Schombeit||2014| ) [hereafter BdJOl, BOS, GB09, Z09, Til, IP13, and 
MS 14 respectively]. These calibrations differ largely by virtue of their underlying SPS treatments, 
including the effects of dust attenuation, bursts of star formation and chemical evolution. These 
effects can also be modelled via other methods for measuring M^/L (e.g. SED fitting), though 
usually at the cost of increased computational expense. The MLCR method suffers no such penalty 
for added model complexity though, making it an expedient, and thus practical (if reliable), route 
to M*/L in most circumstances. 


Despite their widespread use, an exhaustive analysis of the accuracy of MLCRs is still lacking 
(GB09 being the most in depth thus far). BdJOl showed that uncertainties in their adopted treat¬ 
ments of SPS, galaxy evolution, and dust translated into a wavelength-dependent error of 0.1-0.2 
dex in M*/Lwhile large bursts of recent star formation (10% by mass) inflate this range to 0.3-0.5 
dex. Although B03 echoed these conclusions, the BdJOl MLCRs targetted spiral galaxies alone 
(as did MS 14) and the B03 MLCRs excluded treatments for reddening and bursts of star forma¬ 
tion. Z09 argued that omitting the latter effects results in systematically larger M^/L values at 
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blue colours and smaller values at red eolours. Using moek galaxies, GB09 studied a multitude 
of effeets on the reeovery of M*/L, finding, amongst other things, that with suffieient signal-to- 
noise {S/N; >30), individual optical colours performed as well as spectroseopic indices in the 
median, except for galaxies described intermediate-age populations and a bursty star formation 
history (SFH). 


To be sure, effeets other than our assumptions about the SFHs and dust contents of nearby 
galaxies eontribute to the uneertainty budget underlying MLCRs. Z09 eompared MLCRs derived 
from the Bruzual & Chariot ( 2003| hereafter BC03) SPS model and the unpublished 2007 up¬ 
date thereof (hereafter CB07), where the greatest differenee between the two presumably lay in 
their preseriptions for TP-AGB stars. They demonstrated that this phase of stellar evolution alone 
ean alter M^/L predietions by up to 0.1 and 0.4 dex at optieal and near-infrared (NIR) wave¬ 
lengths, respeetively. Til partly explored the effeet of constraining MLCRs with SED fits to real 
galaxies. While their empirical MLCR differs markedly in shape from the one defined by their 
prior-weighted SPS model library (see their Fig. 13), the two appear to be statistieally consistent 
everywhere exeept in the interval g-i = 0.7-1.05. MS14 presented eorreetions to several sets of 
published MLCRs (including B03, Z09, and IP13) for inconsistencies in predicted M^/L through¬ 
out the optieal-NIR wavelength range, presumably tied to uneertainties in stellar evolution, al¬ 
though the eorreetions were solely based on a limited sample of nearby spiral galaxies and the 
degree of ineonsisteney was not quantitatively established. Finally, the ehoiee of IMF is a erueial 
parameter of MLCRs as it eontrols their normalisations (zero-points). For example, M* estimates 
based on a Chabrier or Salpeter IMF differ by 0.3 dex, with the latter being higher. 


Most, if not all, of the systematies inherent to the MLCRs have also been studied in the context 
of SED fitting, all yielding similar results. While much effort has been spent on understanding the 
latter (see |Conroy|2013| and |Courteau et al.|2014| for reeent reviews), the computational expense as¬ 
sociated with even modest parameter spaees severely complieates the investigation of eertain issues 
via SED fitting. Eor instanee, masses derived from integrated photometry are biased to the brightest 
regions within galaxies (i.e. light-weighted) and therefore suseeptible to erroneous M^/L determi¬ 
nations. Two elear eauses of this are interstellar dust or bright stars, whieh ean shroud or outshine 
a potentially large fraetion of a galaxy’s stellar mass. This bias ean be overeome by adding up 
the mass in galaxies pixel-by-pixel. Given the many pixels that loeal galaxies fill with modern 
imaging cameras, measuring M* pixel-by-pixel through SED fits is potentially expensive. Thus in 
their study of spatially-resolved stellar masses, Z09 drew on the MECR method. They found that 
the pixel-by-pixel approach yields masses that are larger than their integrated counterparts by up 
to 0.2 dex, depending on the filters and SPS models used. The Z09 sample however was small 
(9 galaxies) and skewed to late-type spirals (5/9 had Sbe/Se morphologies). An expanded sample 












-4- 


would cement their findings on firmer ground^ 


We wish to address several open issues pertaining to MLCRs for the measurement of M*. 
While earlier studies have focussed on comparisons of theoretical MLCRs against one another 
for various combinations of colours and M^/Lx (the mass-to-light ratio in filter X), we also take 
the approach of comparing them against real galaxies. To this end, we conduct several tests of 
MLCRs using a representative sample of local galaxies drawn from the “SHIVir survey” of Virgo 
Cluster galaxies, for which deep luminosity profiles at optical and NIR wavelengths are available 


(McDonald et al. 2011). We also consider multiple colours at once in fitting for M^/L, whereas 


previous analyses used two colours at most (Z09). We specifically address M* differences due to: 
(i) MLCRs versus SED fitting; (ii) different MLCR sets; (iii) limited SED constraints; and (iv) 
integrated versus spatially-resolved photometric data. 


Our paper is organized as follows. In Section]^ we introduce the SHIVir survey upon which 
our stellar masses are derived. In Section]^ we construct our own MECRs, including a discussion 
of our assumed priors, and test the fidelity of SED-based masses using mock galaxies, which allows 
us to rank the systematics inherent to the MECR method. The latter are investigated and quantified 
(using real galaxies) in Section Einally, our conclusions and suggestions for future work are 
reported in Section 


2. Data 


Thanks to its proximity and richness, the Virgo Cluster enables the simultaneous analysis 
of the stellar populations in all types of present-day galaxies on both spatially-resolved and inte¬ 
grated scales. The SHIVir survey ("Spectroscopy and //-band Imaging of Virgo cluster galaxies"; 
McDonald et al.|20TT | seizes this opportunity by combining optical and NIR imaging for a volume- 
limited sample of Virgo galaxies with apparent R-band magnitudes >16 and located <1.7 Mpc 
(in projection) from the core of the sub-cluster A, taken to be M87. Eike the cluster itself, the 
sample includes the full range of morphological types and is thus representative of the low-redshift 
galaxy population. Eor these reasons, we draw on the SHIVir survey for our tests of the MECR 
method. 


The SHIVir galaxies were imaged in four optical filters {griz) and the NIR //-band. The opti¬ 
cal imaging was downloaded from the seventh data release of the Sloan Digital Sky Survey (SDSS; 


'During the final review of this paper, the manuscript by (Sorba & Sawicki 


2015 


hereafter SS15) expanded upon 
Z09’s findings by measuring pixel-by-pixel masses for 67 nearby galaxies comprising 9 E/SO’s, 4 Irr’s and the rest 
spirals. SS15 stated that the discrepancy between resolved and integrated masses grows with higher specific star 
formation rates and disappears for spatial resolutions >3 kpc. 
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Abazajian et al.||2009). Conversely, the NIR imaging is a heterogeneous colleetion obtained from 


either the GOLDMine database (Gavazzi et al. 2003), Two Micron All Sky Survey (2MASS; Skrut- 


skie et al.]|2006 ), or an extensive observational campaign conducted by our team with Mauna Kea 


telescopes and detectors. SHIVir data products include luminosity profiles, integrated magnitudes, 
effective radii and surface brightnesses, and concentrations for 286 Virgo cluster galaxies in the 
grizH filters. Further details about SHIVir photometry and its associated data products can be 
found in McDonald et al. (2009 201l| and at this website^ 


3. Modelling the Stellar Mass-to-Light Ratio 


We present our adopted method for modelling the stellar mass-to-light ratios of SHIVir galax¬ 


ies (Section 3.2). Recall that a major aim of this work is to assess whether M* derived from 
integrated photometry are consistent with those based on spatially-resolved data. Since the latter 
typically involve modelling M^/L for thousands of pixels per galaxy, our method invariably draws 
upon MLCRs for their expediency. Although several calibrations of MLCRs already exist in the 
literature (B03, Z09, T11, IP13), our experiment calls for the development of our own MLCRs, as 
argued below (Section |3^ . We also develop in Section [3T| mock galaxy tests for the recovery of 
M* based on SED fitting. Finally, in Section [T^ we outline our approach to measuring M* through 
ID light profiles and 2D maps of the stellar mass density (S*) distribution in SHIVir galaxies. 


3.1. Mock Galaxy Tests of SED Fitting 

This work focusses, in part, on modelling M* for a representative sample of local galaxies 
whose light (and thus mass) can be measured in three ways that we refer to as OD, ID, and 2D. 
‘OD’ refers to the end point of a galaxy curve of growth, i.e. the total luminosity of the galaxy 
up to a certain radius. ‘ID’ refers to the profile of galaxy surface brightness as a function of 
galactocentric radius; the brightness at each radius represents the median flux divided by the area 
subtended by the associated isophote. Lastly, ‘2D’ stands for the 2-dimensional image of a galaxy 
whose pixels reveal information about the local stellar population. 

The principal unknowns governing A/*/L for any stellar population are its IMF, age, metallic- 
ity (Z), and reddening. Given the wide range of scales explored in this work, our modelling must 
therefore include populations covering a variety of evolutionary states. Specifically, we require 
a library of theoretical stellar populations spanning extensive ranges in star formation rates and 


^See also http : //www . astro . queensu . ca/virgo 
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histories, chemistries, and dust contents. 


The MAGPHYS library (Multi-wavelength Analysis of Galaxy PHYSical Properties; da 


Cunha et al. 20081 is well-suited for our purposes. This library consists of 50,000 randomly- 


generated stellar and dust emission SEDs that can be self-consistently linked to interpret UV 
through sub-mm imaging of galaxies in terms of their stellar population and ISM properties. We 
disregard the dust emission SEDs since constraining ISM properties is not our current concern. 
Assuming that dust emits at wavelengths >2.5 /im, we use only UV-NIR SED information from 
MAGPHYS. This wavelength range comfortably captures the baseline spanned by the SHIVir 
database. 


The MAGPHYS stellar SEDs can be computed based on either the BC03 or CB07 SPS mod¬ 
els. These SEDs are computed according to SEHs prescribed by five parameters. These are the 
start time (1) and exponential decay rate (2) of continuous star formation, as well as the ages (3), 
durations (4), and mass fractions (5) of superimposed bursts, if any. No prior constraints are made 
as to the number of bursts that underlie any given SED other than that half of them reflect a burst 
that occurred within the past 2 Gyr. 


Attenuation is folded into the stellar SEDs following the Chariot & Eall (2000) two-component 
model for dust in molecular clouds and the ambient ISM. The amount of attenuation suffered by 
each SED is described by two parameters: the total optical depth towards stars <0.01 Gyr-old and 
the fractional contribution of ambient dust to that optical depth (/i). The SED metallicities are also 
allowed to cover the range Z = 0.0004-0.04.(i.e. 0.02-2 Z©). Altogether, eight parameters are used 
to describe the stellar SEDs. To compute each SED, the MAGPHYS software randomly draws val¬ 
ues of these parameters from prior distributions. See |da Cunha et aL] ( |2008| ) for more information 
about these priors. 


The extensive parameter space covered by the MAGPHYS library encapsulates the variety of 
stellar populations within present-day galaxies. Accordingly, this library may be used to evaluate 
the accuracy of M* measurements from SED fitting. Specifically, we wish to know the effects that 


S/N and SED sampling have on these measurements. Using mock datasets from MAGPHYS, da 


Cunha et al.| (2008) were able to recover input stellar masses with very high fidelity in the case 


of an exquisitely sampled SED (i.e. UV through sub-mm fluxes) and a fixed uncertainty of 10% 
in each filter. The authors also claimed that the masses measured from the integrated SEDs of 
three galaxies (NGC 0337, 3521; Mrk 33) remained relatively robust even as UV, Balmer, far-IR 
and sub-mm fluxes were iteratively excluded from their fits. These results, however, were not 
quantified. Moreover, the authors did not quantify the variations of their M* estimates if their 
only available data consisted of optical plus NIR fluxes, a situation that we face here and often 
encountered in the literature. 
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To assess the quality of stellar masses that can be achieved with optical and NIR photometry, 
we have conducted our own mock galaxy tests for various combinations of S/N and SED sam¬ 
pling. For each test, 500 randomly selected MAGPHYS stellar SEDs had Gaussian noise of a fixed 
percentage applied to their fluxes from a select filter suite. Stellar masses were then measured from 
these perturbed SEDs using the Bayesian fitting software included in the MAGPHYS package. 

Figure[2demonstrates the accuracy to which MAGPHYS can recover M* for galaxies and how 
that accuracy might vary with galaxy parameters. The parameters studied include mass-weighted 
age, metallicity, ISM optical depth and age of the most recent burst of star formation (tburst)- Results 
are shown for four cases: (i) SEDs perturbed by 10% noise and sampled with the GAEEX FUV 
and NUV, SDSS ugriz, and 2MASS H filters [blue points]; (ii) as in (i) but with 5% noise [red]; 
and (iii) as in (ii) but with SED sampling restricted to the grizH filters [green]; or (iv) using gi 
filters only [orange]. Cases (i) and (ii) in Fig. [^reflect the situation of a well-sampled stellar SED 
while cases (iii) and (iv) match the data available from the SHIVir survey. Case (iv) examines the 
impact of removing r-, z-, and //-band data, on account of the nebular contamination, low S/N, 
and/or poor sky subtraction that sometimes compromises galaxy photometry at those wavelengths. 

In all four cases, the masses recovered by SED fitting are described by a residual distribution 
having zero mean and dispersion of 0.09-0.13 dex, i.e. random errors of 23-30%. The robustness 
of the residuals in Fig. is remarkable given the wide range in data richness explored. We expect 
that the constraints from well-sampled SEDs should better differentiate themselves from those of 
coarser SEDs as the S/N increases beyond what we have explored. Furthermore, accuracies of 
~0.1 dex can be achieved with the SHIVir filters and S/N = 20. Such a threshold S/N translates 
into flux and colour errors of ~0.05 mag and ~0.10 mag, respectively. The latter error will guide 
the construction of our own MECRs, described in the following section. 

Fig.[^also hints at a possible, but complex, trend in accuracy with burst age, whereby masses 
are maximally overestimated around tburst ~ 1-2 Gyr and underestimated at the youngest and 
oldest tburst (exclusive of tburst )• This trend is consistent with the results of GB09, who found 
that, for an error of 0.1 mag in g-i, the M*/L of mock galaxies described by bursty SFHs and 
young/intermediate-age populations have predicted masses with median residuals of 0.0/0.15 dex. 
It is worth noting that the trend appears weaker for the red points, which have higher S/N. These 
results suggest that bursts are poorly constrained with broadband photometry of marginal quality 
{S/N < 10), no matter how many filters one has data for. Further work aimed at addressing the 
combination of filters (including narrow bands) and S/N required to unravel the contribution of 
bursts to the recent SFHs of galaxies would be useful. 
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Fig. 1.— Recovery of the known stellar masses (M*) of 500 mock galaxies via Bayesian SED 
modeling = fitted M*; = true M*], plotted as a function of mass-weighted age, total 

metallicity (Z/Zq), ISM optical depth (tism) and age of the latest burst of star formation (tburst)- The 


galaxies were drawn from and fitted with the MAGPHYS library, based on the Bruzual & Chariot 


( [2003[ hereafter BC03) stellar population synthesis (SPS) model. The experiment is carried out 
for four combinations of SED sampling and signal-to-noise ratio (S/N): (a) UV to near-infrared 
coverage and S/N = 10 (blue points); (b) as in (a) but S/N = 20 (red); (c) SDSS griz + 2MASS 
H, S/N = 20 (green); and (d) SDSS gi (orange), S/N = 20. The filter combinations examined in 
cases (c-d) focus on the data available from the SHIVir survey ( McDonald et al.||2011 ). Datasets 
(a-c) have been offset vertically for display purposes (denoted by the solid lines). Galaxies which 
have experienced a burst during the most recent timestep (i.e. tburst = 0.0) have been assigned to 
log tburst = 6.3, also for display purposes. The residual distribution is similar in all cases, having a 
mean of roughly zero and dispersion of ~0.1 dex. 
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3.2. Estimating M^/L via Colour Information 

In Section [37T| we quantified the internal accuracy of M* measurements from SED fitting. This 
provides us with a benchmark to evaluate the relative success of our adopted method (MLCRs) for 
measuring M^/L (Section [43] ) . As stated earlier, the appeal of MLCRs lies in their minimal data 
requirements and expediency in all applications. 

Since the pioneering work of BdJOl and BOS, MLCRs have been a popular tool for measuring 
M^/L. Similar to the strategies of Z09, GB09, and MS 14, we do not restrict our modelling to a 
single colour. The motivation for this stems from the fact that every colour has a unique pattern of 
sensitivities to the parameters that determine M*/L. As a result, the MLCR for any colour contains 
scatter (Til), which should be accounted for by modelling multiple colours at once. However, 
using all of the N{N-1)I2 colours that may be measured from an SED sampled with N filters would 
implicitly weigh certain parts of the SED more than others. Moreover, modelling all possible 
colours at once is tantamount to overfitting the data. Eor any given SED then, we only consider 
as many colours as are necessary, when paired with the flux in a single band, to reconstruct the 
original SEE|^ We select those colours whose MLCRs have the least amount of scatter; that is, 
those which are in principle the strongest discriminants of M^/L. The fact that we can model 
several colours at once makes our approach unique from those of Z09, GB09, and MS 14, each of 
which considered two colours at most. 

Eigure [^highlights our multi-MLCR method for measuring M^/L . We show with black lines 
theoretical trends in median M^/Li versus all ten colours that can be measured from the full SHIVir 
SED. These trends were determined from the BC03 version of the MAGPHYS library. The grey 
lines trace the 16* and 84* percentiles of the model distribution for each MLCR, which we take 
to be the l-a uncertainty in M^/Li. The red lines indicate the M*/L,values that we would infer 
for a random MAGPHYS model according to its integrated grizH SED. The MLCRs for g-r, 
g-i, g-z, and g-H provide the tightest constraints on M^/Li for an SED sampled with those 
filters, and so are the ones that we would use in this case; hence the solid red lines in those panels 
while all the others are dotted. Comparisons of the values inferred for SHIVir galaxies 

based on these colours confirms that, on average, they are mutually consistent within la. Thus by 
comparing a set of measured colours - be it for an individual pixel, isophote, or whole galaxy - to 
such MLCRs as displayed in Eig. [^ we can deduce a corresponding set of M*/L, measurements. 
These in turn may be combined to obtain a single A/*/L, value for the object in question. We 
use a weighted average, based on the 16-84* percentile range (linear scale), in order to reduce 


^The fact that a flux is required to turn M»/L into means that the number of colours needed to reconstruct a 
given SED is minimized. For instance, a griH SED is fully described by the i-band flux plus the colours g-r, g-i, 
and g-H. 
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the contribution from colours that are poorer M^/Li discriminants. This is how we use multiple 
MLCRs simultaneously to measure M*/L. 

A potential pitfall of our method is that, in trying to minimize the random error, we increase 
the systematic error above the level ineurred were we to use MLCRs with larger scatter instead. 
To show that this is in fact not the case, we plot in Figure the mean offset of the MLCR-based 
M^/Lj values predicted from each of the 10 eolours accessed by the SHIVir dataset as a function 
of the mean scatter in the assoeiated MAGPHYS/BC03 MLCRs. The offsets are computed with 
respect to the M*/L, inferred from fitting the grizH SEDs with the same model library. The figure 
eonvincingly demonstrates that, despite a couple of outliers, there is a clear eorrelation between 
the two variables. In other words, the precision of a given MLCR is, on average, directly related 
to the accuracy of its predictions and our method is physically sound. We have also investigated 
the impaet of folding systematic errors though our mock galaxy tests (Fig. into our weighting 
scheme and found them to be minimal. 


To apply our method to SHIVir data, we require MLCRs east in terms of both the SDSS and 
2MASS/// bands. Existing MLCRs prove inadequate for our purposes because they either do not 
treat both optical and NIR colours (e.g. Z09) or they are based on simplistic priors (e.g. BOS). We 
also wish to measure M* through E* profiles and maps, in which case our MLCRs must incorporate 
SFHs that are generally not smooth. Finally, our desire to eneompass all morphologies, as in 
SHIVir, favors MLCRs that span a wide variety of stellar populations. For these reasons, we 
derived our own MLCRs, based on the MAGPHYS stellar SEDs. 


We calibrate our MLCRs in terms of M^/Li and all possible colour combinations formed 
from the grizH filters. Although NIR wavelengths are often viewed as yielding the most robust 
M^/L values (e.g. BdJOl), the current uncertainties in NIR fluxes predieted for TP-AGB stars 
make us favour an optical band instead ( Maraston|2005 ; Bruzual||2007[ ). For each relevant colour, 
we calibrate the associated MLCR with the MAGPHYS/BC03 library in the following way. First, 
A/*/L, and the colour of interest are eomputed for eaeh SED in the library, where the former is 
based on the present-day M*. Those SEDs are then grouped by this colour into 0.1 mag-wide bins. 
The median plus the 16* and 84* percentiles of the M*/L, distribution in each bin are then caleu- 
lated, where the pereentiles are used to approximate the Ll-a uneertainties. This straightforward 
procedure yields the desired MECR. Note that the choice of bin width was made based on Eig. 
which shows that optical-to-NIR SEDs having S/N = 20 yield stellar masses of comparable ac- 
curaey to those with better sampling. Eurthermore, surface brightness measurements with S/N > 
20 may be aehieved in the outskirts of galaxies through, e.g., isophotal fitting. Eor colours having 
errors >0.1 mag, we adjust the l-a uneertainties on their M*/L, measurements to refleet the 16* 
and 84* percentiles of the bins in whieh their lower and upper limits fall, respectively. Although 
our approach to such cases favours larger uneertainties, a more aecurate treatment would introduce 
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Fig. 2.— Demonstration of the simultaneous use of multiple eolours to eonstrain the stellar mass- 
to-light ratio (A/*/L) via A/*/L-colour relations (MLCRs; eonstrueted using the MAGPHYS/BC03 
library). The black line in each panel traces the median M*/L, as a function of the particular colour 
plotted along the v-axis, while the gray lines indicate the 16* and 84* percentiles. Here we fo¬ 
cus on the colours that can be measured for galaxies with the SHIVir dataset. Red lines mark 
the M^/Li values that would be measured for a hypothetical population, given its integrated grizH 
SED. Since four colours plus the /-band flux is the minimum amount of information required to 
uniquely describe such an SED, we differentiate the red lines into solid/dotted types according to 
those colours that provide the strongest/weakest constraints on M*/L,. Our procedure then com¬ 
bines the strongest estimates via a weighted average to obtain the overall M*/L, for the population 
in question. While this demo is based on integrated data, it is straightforward to generalize our 
procedure to work with colour proflles/maps instead. 
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Fig. 3.— Absolute mean offset of M^/Li predieted by MLCRs versus the mean seatter in the 
relations. Eaeh data point eorrespnds to a eolour that ean be aeeessed by the SHIVir survey (e.g. g- 
i) and the value along the vertieal axis represents an average over the 280 galaxies eontained 
therein. Offsets were measured with respeet to the M^/Li determined by fitting the grizH SEDs of 
the SHIVir galaxies direetly. Aside from a eouple outliers, there is a elear eorrelation between the 
two variables, implying that the preeision of a given MECR is related to its aeeuraey. 
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precisely the type of eomputational intensity we wish to avoid with our method. 


3.3. Stellar Mass Density Profiles and Maps 


Fig. [^demonstrates how MLCRs reduce the measurement of M* to a trivial algorithm. This 
defining aspect of the MLCR method makes it highly desirable when considering, say, SED mod¬ 
elling for hundreds of thousands of resolution elements per galaxy. Recall that the motivation 
to compute M* from 2D maps stems from Z09’s finding (corroborated by SSI5) that the mass 
assessment is more accurate that way. 


We now describe the additional preparations to infer M* for galaxies based on either multi¬ 
band luminosity profiles or full images. Note that we shall refer to masses derived from integrated 
luminosities, surface brightness profiles, and images as Ml^, and Ml^, respectively. 


For the purpose of comparing these various sets of masses, it is crucial that they be measured 
at a common galactocentric radius. We choose the f-band 23.5 mag arcsec“^ isophotal radius 
(^ 23 . 5 ,i) since it and the enclosed light (m 23 . 5 ,,) yield the tightest scaling relations ( Hall et aL||2012 l 
and are good proxys for the total sizes and luminosities of galaxies of intermediate-to-high surface 
brightnesses. This is demonstrated in the top panel of Figure Q which plots m 23 , 5 j against the 
extrapolated magnitude in the same band {m,ot {) for the SHIVir galaxies. The former were extracted 


from the SHIVir surface brightness profiles, while the latter were measured by McDonald et al. 


( |201 1 ). The points have been colour-coded by galaxy concentration and the solid black lines mark 
equality. Three-quarters of the SHIVir sample have (m 23 , 5 ,, -m,o,,) < 0.3 mag, defined by high and 


intermediate surface brightness galaxies, while departures from m 23 , 5 ,. 


m, 


totj 


grow towards lower 


surface brightnesses, as shown in the bottom panel of Fig. [^ Since our main concern in this work 
is whether multiple colour measurements can be used to reliably constrain M^/L, these departures 
should not affect any of our conclusions. 


3.3.1. 

Since our methods for measuring and Ml^ of a galaxy are intertwined, we shall discuss 
both herein. Both approaches begin with the grizH light profiles from SHIVir. Depending on the 
bands being modelled, we interpolate the profiles (locally) for each galaxy to the largest pixel scale 
therein. No interpolation is required when modelling optical SEDs alone, since the original data 
come from a single source (SDSS; pixel scale = 0.396")- On the other hand, the heterogeneous 
nature of the NIR imaging - having pixel scales in the range = 0.25-1.61"- imply that either 
those or the optical profiles must be interpolated to a common scale when optical-NIR SEDs are 







- 14- 



20 22 24 
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Fig. 4.— (top) Comparison of /-band aperture and extrapolated magnitudes ( 11123.54 nitotj) 
for SHIVir galaxies, where the former have been measured at the 23.5 mag aresee“^ isophotal 
radius. The points are eolour-eoded aeeording to galaxy eoneentration, C 28 [C 28 < 2 .86 (blue), 
2.86 > C 28 < 3.31 (green), 3.31 > C 28 < 3.99 (orange), C 28 > 3.99 (red)] and the solid line marks 
the 1:1 relation, (bottom) The residuals between these quantities (A = 11123 . 54 - 1111014 ) are plotted as 
a funetion of /-band surfaee brightness. Not surprisingly, m 23 . 5 ,, breaks down as a proxy for total 
light at fie 4 > 23.5 mag aresee”^; our eonelusions are not affeeted by this eaveat though. 
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considered. Interpolated quantities inelude the surfaee brightness, elliptieity, and flux enelosed 
within eaeh elliptieal isophote. We then eorreet for Galaetie reddening using the updated version 
of the |Schlegel et al.| ( [1998[ ) prescription by |Schlafly & Finkbeiner| ( [20TT] ). 


In order to match the resolution of our MLCRs (0.1 mag; Section 3.2), we use only light 
profiles with S/N > 20 per isophote, per filter. To meet this requirement, we bin the isophotes for 
all of the profiles of a given galaxy, marching outwards from the center. The bin sizes are increased 
until the S/N threshold is met in all bands. This procedure requires as input the S/N as a function 
of galactocentric radius for each profile. To conform with the adaptive smoothing technique that 


we apply in our 2D approach (Section 3.3.2), we estimate the S/N of each isophote assuming 


background-dominated noise. The radius and surface brightness of each bin are computed as 
weighted averages of the respective input values, where the weights are set to the errors in the 
original light profiles. Conversely, the enclosed brightness of each bin is simply set to the value 
for the outermost isophote contained therein. The differential and integrated S/N for each bin 
are used, respectively, to define the statistical uncertainty in its surface and enclosed brightness. 
Lastly, the size of the central-most bin for each galaxy is set to match the maximum seeing disk 
amongst all bands being modelled. This ensures that the colours in these regions are accurately 


measured. The seeing values for the SHIVir data were taken from McDonald et al. (2011). Note 


that the binning algorithm employed here is an adaptation of the one described in Roediger et al. 
( |2011aD . 


For each SHIVir galaxy, the above preparations yield refined luminosity profiles and non- 
parametric curves of growth for any desired combination of the grizH filters. From the former, we 
measure the necessary colours for each radial bin and compare these data against a set of MLCRs 
to obtain a profile of . Combining this profile with its /-band curve of growth gives us the 
cumulative M* profile and thus the ID mass within the R 23 . 5 .; aperture. At the same time, the 
integrated colours of the galaxy are extracted from the curves of growth for the same aperture. 
These colours are then modelled with the same set of MLCRs to arrive at the OD mass. The major 
distinction between our OD and ID masses is that the latter account for the luminosity and colour 
gradients in our galaxies, whereas the former simply represent their luminosity-weighted sum. 


3.3.2. 2D masses 

The starting point of our 2D approach to M* measurements is the grizH imaging for a galaxy. 
Recall that, while SHIVir’s optical imaging comes entirely from the SDSS, the NIR imaging orig¬ 
inate from a number of sources. For the latter, we use preferentially the ULBCam images for 
several reasons: (i) the field of view (~ 8.5' x 8.5') permits S* to be measured to large radii, with 
good sky sampling; (ii) the photometric depth is high (/in < 24 mag arcsec“^); (iii) the sensitivity 
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and spatial distortions have been mapped to high preeision; (iv) the sub-sample spans all major 
galaxy types; (v) the photometric calibration is homogeneous and robust; and, (vi) the small pixel 
scale (0.25") ensures high spatial resolution for our mass maps. The major disadvantage of the UL- 
BCam imaging is that it lacks an astrometric calibration, which could compromise the accuracy 
of the measured pixel-by-pixel colours. To overcome this, we computed astrometric solution^for 
our optical and NIR images based on the positions of foreground stars. We have used our optical 
images in this operation for improved homogeneity. 


Other steps are required before a galaxy’s colour and stellar mass maps can be computed. 
The first is to resample the selected images to the largest pixel scale therein (if necessaryLike 
Z09, we do not correct the images to a common PSF since this would corrupt their noise properties 
and thus invalidate our later use of smoothing. Next, we mask foreground stars and background 
galaxies that fall within the vicinity of the target. While we generally base our masking upon 
a single image, we make improvements as necessary for light leaked out in other bands. After 
masking, we place five equal-sized boxes (N = 100 pixels) outside the galaxy, in sourceless regions, 
to measure the background and its variation in each image. The background, which is computed as 
the mean pixel value amongst the boxes, is subsequently subtracted off. These steps leave us with 
a set of co-registered, background-subtracted images with uniform resolution. 


As in our ID masses, we prefer that each resolution element included in our 2D masses meet 
a certain S/N threshold. To this end, we use the technique of adaptive smoothing to boost the S/N 
of pixels that fall below this threshold. The appeal of this technique is that it preserves the spatial 
sampling of the original images. However, it comes at the cost of correlating the noise between 
adjacent pixels in the output. Adaptive spatial binning techniques, such as Voronoi tesselation 
( |Cappellari & Copm||2003| ), avoid this problem but produces images with irregular renditions of 
the 2D information. Since these irregularities would complicate our modelling of the pixel-to-pixel 
colours and extraction of 2D masses (described below), we prefer to raise the 5/A of pixels in our 
images, where necessary, via smoothing. 


We smooth our images with the ADAPTSMOOTH code ( Zibetti|2009 l, assuming background- 
dominated noise. While this assumption ignores the Poisson contribution from the source, the 
pixels that require smoothing are generally located in regions of low surface brightness, where 
the background is still dominant (S. Zibetti, priv. comm.). Moreover, determining the source’s 
contribution to the flux in each pixel requires knowledge of the detector gain, which is often poorly 


"^http://nova. astrometry.net/ 

^This is done with the SWarp software package - http;//www.astromatic.net/software/swarp The output from 
SWarp has fixed dimensions for each band and is centered on the galaxy, so that adequate sky samples may be later 
chosen and pixels are properly registered between images. 
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constrained for stacked images, as in our SHIVir NIR data. We set the noise for each image as the 
rms variation in the background level and threshold S/N of 20, as per the resolution of our MLCRs. 
ADAPTSMOOTH is initially run on selected images to obtain a set of smoothing masks that indicate 
the size of aperture each pixel requires in order to meet the desired S/N threshold in all bands. 
These masks are then maximally stacked to generate a master mask. We re-run ADAPTSMOOTH 
on the original images with this master mask in order to apply a consistent smoothing solution to 
every pixel across all filters. Since smoothing corrupts the noise properties of the output images, 
we conservatively assume that the 5/A^ in each pixel is the threshold value, thus maximizing the 
errors. For all ADAPTSMOOTH runs, we limit the smoothing aperture to a radius <20 pixels. Any 
pixel requiring a larger aperture in any band is clipped within all output images. 


After smoothing, the images are photometrically calibrated at each pixel. Calibration zero- 
points to the AB magnitude system were kindly provided by M. McDonald. We also correct the 
images for Galactic extinction, according to Schlafly & Finkbeiner ( 2011| ). These final prepara¬ 
tions leave us with a set of masked, background-subtracted, and calibrated galaxy images which 
have a common spatial sampling, registration, and minimum S/N per pixel. 


With calibrated surface brightness maps in hand for a select set of filters, we construct the 
stellar mass map of a galaxy by first computing the minimum number of colours for each pixel 
(Sectionand transforming those data into a map of M*/L, through a set of MLCRs. The mass 
map follows by simply multiplying the f-band and maps together. A 2D M* is then extracted 
from the mass map by a curve-of-growth analysis. Specifically, we integrate the mass map using 


a fixed set of elliptical apertures (performed using the XVISTA software package; see Courteau 


1996[ [McDonald et al.||2011[ [Hall et al7||2012| ) that match the ellipticity and position angle of the 


aperture from which the OD/ID masses were measured. This ensures that extraction methods do 
not introduce any observed differences between the three versions of our M* measurements. 


4. Results & Discussion 

We have described in the previous section our method for modelling the colours of galaxies 
in terms of their underlying mass-to-light ratio. We can now use this method to quantify the 
systematics in M* measurements obtained via MLCRs. While previous studies have addressed 
various aspects of this topic from a theoretical (idealised) perspective, we use a representative 
sample of observed local galaxies to focus on the following sources of scatter: (i) SED fitting 
versus MLCRs; (ii) choice of priors; (iii) choice of constraining data; and (iv) number of spatial 
dimensions considered. Testing MLCRs with real galaxies should reveal systematics, if any, that 
are not typically accounted for by the assumed priors. Note that our discussion below is largely 
based on the BC03 MLCRs. 
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4.1. The Correspondence Between MLCRs and SED Fitting 


In Section 3.1 we used mock galaxies to assess the accuracy of M* measurements from SED 
fitting under various conditions of data sampling and quality. Assuming perfect model priors, we 
found that the scatter in recovered M* is approximately constant, with ~0.1 dex, regardless of the 
filter combination within the UV-NIR range and whether S/N = 10 or 20 is achieved (see Fig. [T]j^ 
Since SED-based M^/L measurements should be less biased than those from MECRs (Til), by 
design, this result provides us with a useful benchmark for evaluating the quality of MECR-based 
results, which is the aim of this section. 


We take two complementary approaches to test the correspondence between SED fitting and 
MECRs, both of which focus on the recovery of M*. The first approach uses mock galaxies. The 
advantages of this approach are that: (i) the mass to be recovered is known in each case, and (ii) 
randomly-drawn mock galaxies should populate a significant fraction of the observational space 
spanned by a given MECR set. Therefore, mock galaxies offer a most direct way to uncover any 
biases introduced by MECRs. 

For our mock galaxy sample, we draw 500 stellar SEDs from the MAGPHYS/BC03 library. 
In Figure we assess the accuracies of the MECR-based masses and trends therein with galaxy 
parameters under the same scenarios of data breadth and quality as assumed for the green and 
orange points in Fig.[^(i.e. grizH and gi band combinations, S/N = 20). For simplicity, we use the 
same layout and colour scheme as before. 

Fig. shows that our MECR method yields residual distributions that closely resemble those 
seen in Fig. and, as before, we see hints of trends in the residuals as a function of metallicity 
and tburst- Moreover, the residuals from our method have nearly identical properties to those that 
were found for SED fitting, with a mean value of (essentially) zero and dispersion of 0.13-0.14 
dex. Similar results are obtained when other combinations of the SHIVir bands are used (i.e. giH 
and griH). This complements the work of GB09, who found that M*/L could be satisfactorily 
recovered for mock galaxies using the g-i colour, except those characterized by intermediate-age 
populations and bursty SFHs (where M^/L was biased to higher values by 0.15 dex, in the me¬ 
dian). For galaxies well-described by the MAGPHYS priors then, we infer that there is essentially 
no difference between the masses measured via SED fitting or MECRs. Furthermore, Fig. [^in¬ 
dicates that minimal observational constraints are needed to achieve an accuracy of 33% in the 
M* measurements for such galaxies. 

The correspondence between SED fitting and MECRs for mock galaxies was to be expected 


®In particular, measured from SEDs sampled with only the gi filters fare equally well as those determined using 
filters throughout the UV-NIR wavelength range. 
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Fig. 5.— As in Fig.[^but for M* measured via MLCRs (rather than SED fitting) and SEDs sampled 
with the grizH (green) and gi (orange) filters, and S/N = 20. In both cases we find that the residuals 
have a distribution with zero mean and 0.13-0.14 dex of dispersion. Similar results are found for 
the giH and griH filter combinations. 
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Fig. 6.— (top) M* estimates from our multi-MLCR method versus those from Bayesian 

SED modeling Both sets of masses are derived using the MAGPHYS/BC03 library and 

the integrated grizH SEDs of SHIVir galaxies, (bottom) Residuals, eomputed as A = log - 

log M* are plotted as a function of g-/ colour. The correspondence between the two techniques 
is clearly excellent over the full range of colours probed by our dataset (A = -0.01 dex; a = 0.05 
dex). 
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since we drew their SEDs from the same SPS model library that was used to build our MLCRs. Our 
second approach to testing this correspondence overcomes this bias by working with real galaxies. 
In all likelihood, real galaxies populate M*/L-colour space differently than our synthetic models, 
in which case the two methods should yield discrepant results, as demonstrated in T11. A corollary 
to this is that a sample of real galaxies may not span the same range of colours as our models, in 
which case our second approach can only inform us of the correspondence between SED fitting 
and MECRs over a limited interval. SHIVir galaxies, for instance, may be biased to redder colours 
since they inhabit a cluster environment. Still, the outcome of this second approach complements 
the previous one and will provide some context to forthcoming tests. 

Eigurej^ compares M* for SHIVir galaxies obtained from MECRs and SED fitting. Note that 
this comparison assumes that SED fitting yields the most accurate masses for real galaxies. Both 
sets of M* were determined from the integrated grizH SEDs of these galaxies - the results for 
cases involving other bandpass combinations are rather similar. The top panel of Eig. [^indicates 
that the two methods are in excellent agreement, with a mean residual of -0.01 dex and 0.05 dex 
of dispersion, while the lower panel shows that the residuals are independent of galaxy colour, 
although the sampling for g-i < 0.6 is admittedly poor. This is at odds with Til, who emphasised 
the bias imparted by MECRs on M*, but our use of multiple colours explains this difference. Eor 
instance, if we used the g-i colour only, we would tend to overestimate M^/Lt, in agreement with 
Til (see Eig.[Tg. 

The above tests demonstrate the two following points about our multi-MECR method: (i) 
like SED fitting, it recovers the known M* of mock galaxies to better than ±0.15 dex; and (ii) it 
reproduces A/* from SED fitting to real galaxies to within ±0.05 dex. The increased accuracy in the 
latter case is likely caused by the smaller colour range spanned by the SHIVir galaxies, compared 
to our mock dataset. These findings indicate that our multi-MECR method may be counted on to 
yield fairly accurate results. 


4.2. The Underlying Priors of MLCRs 


We have shown in the previous section that MECRs produce, on average, equivalent results to 
SED fitting. However, this does not entail that MECRs and how they are used are completely free 
of systematics. Since an MECR set is constructed from an SPS model library, the priors underly¬ 
ing the model ingredients or how the models are implemented can contribute significantly to the 
error budget. Specifically, any SPS model library parameter that affects the relationships between 
colour, age, metallicity, reddening, and M^/L will govern our results. Eor instance, omitting a 
burst component from model SEHs will bias M^/L to higher values. Given the well-studied role 
that model priors play in SED-based M* (Conroy 2013[ Courteau et al. 20141, we now examine 








- 22 - 


how these effeets are manifested in MLCR-based M*. We omit the ehoiee of IMF from the follow 
ing sinee its effeet, while signifieant (A M* ~0.3 dex), is well understood in most eases (e.g. 
Dokkum|2008[|Longhetti & Saracco|2009[ [Courteau et al.|2014[ ). 


van 


A straightforward approach to assess the effect that priors have on MLCR-based M* is to 
compare MLCRs calibrated with different model libraries. This is achieved in Figures |7]j^ where 
the MLCRs are either of our own creation or drawn from the literature. These plots are structured 
similarly to Fig. [^though the first focusses on optical colours accessible through SHIVir and the 
second on optical-NIR colours. Our own MLCRs correspond to calibrations based on the BC03 


and CB07 versions of the MAGPHYS library, as well as the SPS model of Conroy et al. (2009 


hereafter FSPS). FSPS provides the proper machinery to create libraries like MAGPHYS. With 
it, we have computed 50,000 SEDs which adhere as closely possible to the MAGPHYS priors. 
Adhering to these priors allows us to gauge the impact that the choice of SPS model has on MLCRs. 
Two exceptions to our adherence include the shape of the effective absorption curve for birth cloud 
dust and the treatment of bursts of star formation. In MAGPHYS, the optical depths of ambient 
and birth cloud dust are parameterized as power laws in wavelength to different exponents (-0.7 
and -1.3, respectively), while bursts occur randomly throughout each model’s lifetime and can last 
anywhere from 0.03-0.3 Gyr. In FSPS, both dust components are simultaneously described by a 
single power-law relation, albeit with an adjustable exponent, while only one instantaneous burst 
can be added to any model’s SFH. For those readers interested in incorporating our BC03 and FSPS 
MLCRs into their own work, we present tabulations of linear least- squares fits to these relations 
in the Appendix. 


For the literature MLCRs in Fig.|^ we show the calibrations of B03, Z09, Til, and IP13. We 
omit the MLCR sets of BdJOl and MS 14 given their Johnson-Cousins band calibrations, whereas 
our focus is on the SDSS -i- 2MASS suite. GB09 did not tabulate the MLCRs used in their work but 
we suspect that they are well-matched to those of Z09. The number of literature MLCRs shown 
in Fig. |7] depends on the colours examined here. For instance, the Til MLCR was only defined in 
terms of g-i. Whenever possible, we have restricted the range of colours spanned by the literature 
MLCRs to those advocated by their respective authors. This is why, for example, the Til and 
IP 13 MLCRs do not span the full baselines of the panels in which they appear. An important 
distinction about the trends shown in Fig.is that those from B03 and Til were constrained by 
their authors using SED fits to real galaxies, while all others are purely theoretical. Eastly, the 
B03 relations have been reduced by 0.15 dex to correct their choice of a “diet-Salpeter IME" to a 
Chabrierl ( [20031 ) IME. 

In Eig. [^we demonstrated that, in the case of our BC03 MECRs, M*/L, follows a largely 
monotonic relationship with each one of the colours defined by the grizH bands. Eigs. |7][^show 
that this behaviour is likewise reflected, albeit to different degrees, by other sets of MECRs. The 
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Fig. 7.— As in Fig. [^but now showing the MLCRs predicted for SHIVir optical colours from 
libraries of different SPS models (BC03; CB07; [Conroy et al .112009] hereafter FSPS) or related 
work in the literature ( [Bell et aL]|2003[ jZibetti et~aL 2009 [ Taylor et'aL]|2011[ [Into & Portinari 


2013 hereafter B03, Z09, Til, and IP13, respectively). Note that the B03 and Til MLCRs are 


constrained empirically by SED fits to real galaxies while all others are fully theoretical - the 
BC03, CB07, and FSPS relations are of our own making. Also, the B03 relations have been 
reduced by 0.15 dex to change the underlying initial mass function (IMF) to the Kroupa/Chabrier 
form used by all other MLCRs shown here and the IP 13 relations correspond to their exponential 
star formation history model. 
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Fig. 8.— As in Fig.|7]but now comparing MLCRs predicted for SHIVir optical-NIR colours, where 
available. 
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reason for this is simple: as the age, metallieity, and/or reddening of a stellar population inereases, 
so do its M^/L and eolours. 


Fig. [^highlights that the best agreement between the optieal MLCRs oeeurs with g-r and 
g-i, while all other eolours are eharaeterized by notable discrepeancies. Not surprisingly, the 
CB07 and Z09 MLCRs eompare very well in all panels since they were both derived from the 
MAGPHYS/CB07 library. BOS notwithstanding, these two calibrations also predict the lowest 
A/*/L, values for any colour; we therefore focus the following comparisons around the CB07 
MLCRs. All other MLCRs, except those of BOS and Til, usually agree better with the CB07 cal¬ 
ibration towards redder colours (i.e. older, more metal-rich, and/or highly-reddened populations), 
but the improvement for IP IS is admittedly marginal. Towards bluer colours (i.e. younger, more 
metal-poor, and/or dust-free populations), these same MLCRs can differ from the CB07 calibration 
by as much as ~0.S dex. The FSPS MLCR, on the other hand, predicts a higher M^/Lt by 0.5 dex 
for the bluest colours. The absolute discrepancy between the Til and CB07 MLCRs falls within 
the range <0.S dex at all colours. 


The above observations from Fig. [^indicate that the systematic error in MLCR-based M^/L varies 
as a function of colour. The cause of this colour-dependent error must be due, in part, to the dif¬ 
ferent physics encoded by their underlying priors. It is noteworthy that the BOS and T11 MLCRs 
possess the shallowest slopes, since these are the only relations shown in Fig.jTJthat were calibrated 
using real galaxies. To what extent this difference reflects a genuine flaw in theoretical MLCRs is 
unclear given the marked differences in the priors between them and the empirically-constrained 
MLCRs. For instance, T11 did not include bursts in their model SFHs and BOS limited the strength 
of such events to <10% by mass, whereas the MAGPHYS library incorporates a variety of bursty 
SFHs. Furthermore, BOS omitted reddening from their models, while Til and MAGPHYS in¬ 
voked the Calzetti et al. (2000) and Chariot & Fall ( [2000 ) extinction laws, respectively. Teasing 
out the relative contributions of priors and the real galaxy distribution in the A/*/L-colour plane 
to the slopes of the BOS and Til MLCRs is non-trivial and lies beyond the scope of this work. 
Furthermore, this comparison between empirical and theoretical MLCRs can only inform us on 
the systematics of modelling integrated light. Investigations like BOS and Til, though based on 
spatially-resolved data, have yet to be performed. 


Turning to Fig. we find that the systematic errors incurred by using optical-NIR colours 
are often much worse than those associated with optical colours alone. This contradicts the pop¬ 
ular belief that NIR data yield the most accurate M^/L measurements. We confirm Z09’s claim 
that the optical-NIR MLCRs for the MAGPHYS/BCOS library tend to exceed those of the MAG- 
PHYS/CB07 library, often by as much as O.S dex. Since the main change to the CB07 model was 
to increase the NIR flux of TP-AGB stars ( |Bruzual|20071 ), we attribute this difference to the role of 
this stage of stellar evolution to the error budget. Supporting this interpretation is the fact that these 
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SPS models are implemented identieally in the MAGPHYS library, and therefore any differenees 
between their MLCRs must be rooted in their most basie ingredients. 


Although signifieant, the differences between the BC03 and CB07 optical-NIR MLCRs pale 
in comparison to those between the CB07 and FSPS MLCRs, which can exceed 0.6 dex. To help 
identify the source of this discrepancy, we show in Figure the MLCRs corresponding to the 
simple stellar populations (SSPs) of the BC03, CB07, and FSPS models, at solar metallicity. Since 
the CB07 and FSPS models incorporate the same isochrones (Marigo & Girardi|2007[|Marigo et ah 


20081, we expect that the MLCRs for their SSPs should follow each other. Instead, the CB07/SSP 


MLCRs abruptly diverge from the FSPS/SSP MLCRs to redder optical-NIR colours at an age 
of ~0.1 Gyr, and remain well-separated for another ~ 1.5-1.6 Gyr. During that time period, the 
FSPS/SSP MLCRs actually follow the BC03/SSP relations. These behaviours qualitatively persist 
at lower metallicities as well (e.g. Z = 0.004). The above age range is well-known as the period in 
which TP-AGB stars are present in SSPs. Therefore, the TP-AGB treatment must be contributing 
to the observed differences between the CB07 and FSPS optical-NIR MLCRs too. However, Fig.[^ 
shows that the FSPS MLCRs also lie above the BC03 MLCRs, despite the similar trends exhibited 
by their SSPs in Fig. We suggest then that the much burstier nature of the MAGPHYS SFHs 
helps drive the differences between the FSPS and BC03/CB07 MLCRs as well. The flat prior 
on the burst mass fraction in MAGPHYS and the fact that half of its SFHs incorporate a recent 
burst suggest that our BC03/CB07 MLCRs are biased to stellar populations with younger mass- 
weighted ages than our FSPS MLCRs. Our suggestion is consistent with the rank order of these 
MLCRs in Figs.|7]j^ in which the BC03/CB07 MLCRs predict a lower M*/L, at bluer colours than 
do the FSPS MLCRs. Caution must therefore be taken when interpreting M* measurements that 
are based on optical-NIR MLCRs. 


4.3. The Impact of SED Sampling 

Figs. |7]j^ highlighted the systematic penalty incurred with the MLCR method through the 
choice of an SPS model library. Particularly striking was the inflated disagreement between MLCR 
sets when optical-NIR colours are included. However, it is conceivable that, within the context of 
a single MLCR set, the M*/L measured for a galaxy would also vary as a function of the data 
considered. One reason could be the fact that not all colours are equally sensitive to the parameters 
that establish M^/L, such that the use of multiple colours might enhance one’s constraining power. 
Another could be that the fluxes predicted by SPS models within one or more filters are inaccurate, 
which likewise would affect all colours associated therewith. Given these concerns, we address in 
this section whether or not the degree of SED sampling contributes meaningfully to the systematics 
of MLCRs. 
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Fig. 9.— As in Fig. but now comparing the MLCRs based on SSPs of solar metallieity drawn 
from the BC03, CB07, and FSPS SPS models, and for all SHIVir eolours simultaneously. Similar 
trends are also found amongst these SSPs at lower metallieities (e.g. Z = 0.004). 
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To assess the impaet of SED sampling, we shall test the robustness of our multi-MLCR 
method under a variety of realizations of the SHIVir data. We eompare in Figure [T^ sets of 
A/* measured for these galaxies via our BC03 MLCRs and SEDs sampled with the following filter 
combinations: gi, giH, griH, and grizH. The layout of this figure is identical to Fig.|^ and we use 
the same combinations of colours as before in the fits for each SED realization. The gi-, giH-, 
and gn//-based M* are plotted in the top panel as a function of the gr/z//-based A/*, which we 
take to be the most accurate (see Fig. [^. In the bottom panel we plot the residuals between the 
former and latter as a function of the g-i colour. 


Fig. [T^ demonstrates that the grizH-hd&ed masses of SHIVir galaxies are reproduced well by 
those using the griH and giH band combinations instead. The g/-based values, on the other hand, 
are biased to slightly higher values, by 12% on average, implying that the g-f colour systematically 
overestimates A/*/L,, in agreement with the findings of Til. In the bottom panel, we see that the 
residuals from the various realizations of the SHIVir dataset are robust to galaxy colour but, as in 
Fig.[^ this conclusion is tentative for systems with the bluest colours {g-i < 0.6). 

The robustness of A/* for SHIVir galaxies is explained by how we combine A/*/L, values when 
multiple colours are fitted. Recall that we construct a master Af*/L, for each object by computing 
a weighted mean of the values that were assigned to its colours. Our method then naturally assigns 
greater weight to those colours which have less scatter in their MECRs. In the case of Fig. [T^ 
the colour that maintains the robustness amongst the various sets of M^, is g- i. Indeed, Fig. 
demonstrates (via the dashed lines) how the scatter in A#*/L, is minimized for the g-i MFCR, 
relative to the g-r, g-z, and g-H MECRs. 


Since most sets of A/* presented in Fig. 10 are mutually consistent, one may be tempted to 
conclude that the r and z bands do not increase the constraining power of our method. While it is 
true that including these bands does not change A/*/L, significantly, a positive statistical impact is 
made on the assigned errors. The mean error on A^*/L, with the gi, giH, griH, and grizH SEDs 
of SHIVir galaxies, is 0.29, 0.22, 0.18, and 0.15 dex, respectively. Thus while the g-i colour may 
be counted on to reasonably determine Af* when other photometry is scarce and/or unreliable, the 
desirable combination of high accuracy and precision is achieved by applying MECRs to well- 
sampled SEDs. 


4.4. Must Mass Estimates be Spatially-Resolved? 

We now examine the claim that Af* is more accurately measured on the basis of spatially- 
resolved data than integrated fluxes (Z09). According to this claim, the method by which one 
measures a galaxy Af* is most relevant if a significant amount of it (e.g. >10%) is hidden by 
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g-i (mag) 

Fig. 10.— (top) Comparison of M* estimates from our multi-MLCR method for various real¬ 
izations of the SHIVir dataset. Masses derived from the filter eombinations gi, giH and griH are 
plotted against those measured from the grizH bands, (bottom) Residuals, eomputed as A = logM;^ 
- logMf"^^, are plotted as a funetion of g-i eolour. The mean and dispersion of eaeh residual dis¬ 
tribution are quoted in the top panel. These values indieate that our method is modestly robust to 
SED sampling. 
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bright stars and/or dust. Higher-dimensional approaehes to measuring M* presumably overeome 
this bias by assigning eaeh galaxy resolution element a M^/L tailored to its eolours. In other 
words, allowing M^/L to vary spatially within a galaxy may be essential to ensure an aecurate 
determination of its M*. 


The need for a 2D approaeh to measuring M* was first demonstrated by Z09. Using optieal 
and NIR imaging, they found that the of spiral galaxies eould be lower than their by as 
mueh as 0.15 dex. Conversely, the 2D approaeh produeed hardly any ehange in M* for the lone 
early-type galaxy in their sample. The authors thus elaimed that the diserepaney between and 
inereases towards later galaxy types. The Z09 sample, however, was small and biased to late- 
type spirals. We have improved upon this limitation by expanding the work of Z09 to a larger, more 
representative set of loeal (SHIVir) galaxies (see also SS15). Similarly, the use of a ID approaeh 
(via luminosity profiles) may also yield improved M* measurements, by aeeounting for the eolour 


gradients which many galaxies are known to possess Roediger et al. (2011a). A ID approach is 
also more straightforward and less computationally-expensive than its 2D cousin. We therefore 
explore the relation between and in this section as well. 

Figure [TT] compares the OD and ID M* for SHIVir galaxies based on SEDs sampled with 
either the gz, giH, griH, or grizH bands. The correspondence between the OD and ID values is 
clearly excellent in most cases, in the sense that there is no net bias, as evidenced by the residual 
statistics. The major outliers are understood as being due to cases of extreme colour gradients, 
where a positive gradient would lead to > M^^and vice versa. Omitting the outliers, we infer 
that radial variations of M*/L, within galaxies, on average, bias OD measurements by no more than 
±20%. Several explanations could apply here: (i) most galaxies have mild colour gradients; (ii) as 
in (i) but only for those colours which possess steep MLCRs (e.g. g-r), or (iii) M^/Li variations 
occur in regions of galaxies which do not enclose large fractions of their total light. We propose 
that the third explanation applies most often, as earlier work by our group ( Roediger et al. |20 11 a I 
demonstrated that significant colour gradients are commonly found within the SHIVir sample. 


Fig. ^Tj suggests that colour gradients do not affect M* measurements for our galaxies to a 
large degree. In addition to the radial coordinate though, galaxy colours can vary with azimuth and 
asymmetrically (i.e. pixel-by-pixel). ID masses are exempt from these other possible dependencies 
by virtue of the isophotal fitting that is used to measure luminosity profiles of galaxies. 2D masses, 
however, are sensitive to all three. To address the significance of the latter two dependencies, 
and the veracity of Z09’s conclusions, we have derived 2D masses for 18 nearly face-on SHIVir 
galaxies, which is twice the size of Z09’s sampleQSome basic properties of our 2D sub-sample are 


^Computational constraints limited us to this sub-sample. We will present the mass maps for all the ~280 SHIVir 
galaxies in a forthcoming publication. SS15 have also recently computed 2D masses for 67 nearby galaxies. 









Fig. 11.— M* estimates from the integration of stellar mass density profiles versus those 

from integrated photometry Both sets of masses have been derived for various realizations 

of the SHIVir dataset (gi; giH-, griH; grizH). The mean residuals and dispersions (tabulated at 
lower-right) suggest that, for our sample, eolour gradients do not introduee a persistent, signifieant 
bias inM°^. 
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listed in Table The morphologieal distribution (eolumn 5) show that this sub-sample represents 
the present-day galaxy population quite well, spanning the giant early-types, spirals, and both 
flavours of dwarves (gas-rieh and -poor). Their OD M* should thus eneompass the full range of 
biases mentioned above, if indeed present. 

The OD and 2D galaxy M* are eompared in FigureAs for Fig. we have measured 2D 
A/* (see Table[^ for a variety of colour combinations, which amounts to iteratively dropping filters 
from the full SHIVir suite (grizH). Unlike our OD-ID comparison, we plot the ratio log(M^^/M°^) 
instead, where galaxies with lie above the horizontal line that marks equality. Plotting 

this ratio allows us to order galaxies from early-types on the left to late-types on the right, and 
thus to gauge the existence of any correlations with morphology. Note that each vertical cluster 
of blue, green, orange, and red points represents a single galaxy and the dwarf galaxies in our 2D 
sub-sample appear as the four leftmost members of the early-type group and the four members of 
the "Im/BCD" group. Arrows represent cases where only lower limits could be placed on the 2D 
masses, while hollow points denote cases in which OD masses are unreliable. The former occur 
when ADAPTSMOOTH cannot build enough S/N out to i? 23 . 5 .i and the latter when one or more 
colours exceed the ranges of our MLCRs. The error bars shown in the lower-right corner represent 
the average uncertainty in the OD masses for these galaxies, and should thus be regarded as lower 
limits. 


Fig. 12 shows that our 2D approach produces mixed results with respect to OD values, de¬ 
pending on the data involved. most often exceeds when the SEDs of pixels are sampled 
with the gi band combination. In the median however, this excess is quite small for early-types 
(dwarves -i- giants; 0.02 dex) and zero for spirals, the reason being that many members of these 
groups lie about the line. While late-type dwarfs possess amongst the highest values of 

log when the gi bands are used, the small size of their group prevents us from drawing 

firm conclusions. Nevertheless, the offsets seen for some of our Im/BCD’s agrees with the fact that 
these galaxies often host pockets of recent star formation, which would far outshine their underly¬ 
ing (low surface brightness) stellar disks/spheroids. In other words, we expect the integrated light 
of these galaxies to be biased blue and, accordingly, low M^/L. 


Although the offsets between and are typically small for most galaxy types when 
using only the gi bands, they can be most significant on a case-by-case basis. For instance, VCC 
1437, an elliptical galaxy with a very blue center (most likely due to recent star formation), has 
> A/”^ by 0.15 dex, while three of the Im/BCD galaxies in our 2D sub-sample have absolute 
log values of >0.1 dex. Our foremost conclusion regarding A/°^ measurements based 

on optical information alone is that internal variations in A7*/L may be captured by including an 
additional random component in the error budget for all populations. Our best estimate of the size 
of this component is 0.05 dex for early-types and 0.03 dex for spirals; these values can be refined 
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Table 1. Bulk properties of Virgo eluster galaxy sub-sample used for 2D mass map experiment' 


NGC 

ID 

(1) 

vcc 

ID 

(2) 

02000 

(deg) 

(3) 

<^2000 

(deg) 

(4) 

Type^ 

(5) 

(mag) 

(6) 

C 28 

(7) 

He 

(mag arcsec“^) 
(8) 

Re^ 

(kpc) 

(9) 

— 

0389 

185.0138 

14.9614 

dSO 

-18.32 

3.690 

22.16 

1.454 

- 

0751 

186.2012 

18.1950 

dSO 

-17.68 

3.374 

22.00 

0.994 

- 

0788 

186.3200 

11.6053 

dE 

-15.98 

2.992 

23.26 

0.970 

- 

0816 

186.4000 

15.8478 

dE 

-17.70 

2.665 

24.46 

3.210 

- 

1437 

188.1396 

9.1736 

E 

-16.96 

4.055 

21.00 

0.437 

- 

1440 

188.1392 

15.4153 

E 

-17.87 

5.797 

22.12 

0.914 

4503 

1412 

188.0258 

11.1764 

SO 

-20.76 

4.701 

20.83 

2.308 

- 

1358 

187.8467 

17.2064 

Sa 

-16.91 

2.679 

23.23 

1.361 

4584 

1757 

189.5746 

13.1100 

Sa 

-18.65 

3.464 

22.05 

1.744 

4407 

0912 

186.6342 

12.6111 

Sab 

-19.64 

3.411 

22.07 

2.919 

4394 

0857 

186.4812 

18.2142 

Sb 

-21.03 

4.239 

22.01 

4.284 

- 

2023 

191.3833 

13.3325 

Sc 

-17.80 

1.943 

23.33 

2.377 

- 

0510 

185.4738 

15.6458 

Sd 

-17.42 

3.391 

23.55 

1.916 

4519 

1508 

188.3758 

8.6547 

Sd 

-19.50 

2.897 

22.12 

2.775 

- 

0583 

185.6875 

15.5022 

Im 

-16.83 

2.384 

24.04 

1.758 

- 

0980 

186.7967 

15.8967 

Im 

-17.52 

2.488 

23.60 

2.373 

- 

1725 

189.4229 

8.5586 

BCD 

-17.23 

2.052 

23.30 

1.777 

- 

1804 

189.9171 

9.3989 

BCD 

-16.35 

3.427 

23.70 

1.324 


^The parameters M, C 28 , He, and Re correspond to /-band values. 
^Drawn from the NASA Extragalactic Database. 


^Estimated assuming a common distance to all Virgo cluster galaxies of 16.5 Mpc (i.e. m-M 
31.087 mag; Mei et al.|2007 1. 
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Fig. 12.— Comparison of M* values derived from integrating the S* maps for 18 SHIVir galaxies 
(MI^) versus The eomparison is plotted as the logarithmic difference between these two sets 
of masses [log(Ml^/M^^)] as a function of morphology and is carried out for SEDs sampled with 
the filter combinations gi, giH, griH, and grizH. Galaxies are ordered from early-types (dE, dSO, 
E, and SO) on the left, to spirals (Sa-Sm) in the middle, and gas-rich dwarves (Im/BCD) on the 
right, with dashed lines separating each of these groups. The solid line denotes and 

representative error bars are shown in the bottom-right corner. 
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using a larger and more eomplete sample in future work. 


2D masses that inelude NIR information behave similarly to those based on gi bands, exeept 
for two important differenees: we rarely find positive values of log and few galaxies 

eluster about the line. The first differenee eould arise from either an increase in or a 

decrease in (or both) as we model more bands than just g and i. Fig. [^shows that, on average, 
A/°^ measured from the band combinations giH, griH, and grizH are mutually consistent (to within 
±0.03-5 dex), while g/-based masses are ~13% higher. This behaviour of can therefore not 
explain the observed changes in log the onus must instead lie with A#^^. Taking ratios 

of between the datasets considered in Fig. 12 indeed shows that typically decreases by 


~30% when the SEDs of pixels are sampled with either the giH, griH, or grizH filters, compared 
to the case of gi SEDs. Including NIR fluxes in 2D mass measurements thus systematically reduces 


galaxy-wide mean A/*/L, values (as in the OD case; Fig. 101 and mass maps derived for each of the 
band combinations gr, gz, and gH confirm that the lattermost yield the smallest values of M^/Li per 
pixel, on average. The discrepancy between optical and optical-i-NIR masses from our 2D approach 
supports the suggestion of Til that predictions of current SPS models lack consistency across the 
optical-NIR domain. Our own results highlight the need for further tests of the MLCR method (and 
SPS models in general), preferably on resolved populations, where M^/L may be directly inferred 
from colour-magnitude diagrams. 


M 


ID 


Unlike with our gz-based measurements, the lack of galaxies found in the vicinity of the 
= line in Fig. 12 when NIR data are considered implies that OD values are biased higher 


than their 2D counterparts in such cases. For early-types and spirals, this bias remains robust 
amongst the giH, griH, and grizH band combinations, falling between -0.05 and -0.07 dex. The 
bias for gas-rich dwarfs changes significantly between these cases, however, from -0.07, to -0.09, 
and -0.15 dex, respectively. In addition, log measurements scatter about these median 

offsets at roughly the 0.04, 0.05-0.06, and 0.10-0.15 dex level for spiral, early-type, and gas-rich 
dwarf galaxy types. Again though, these estimates must be revisited when larger, homogeneous, 
and more complete samples of high-quality optical -i- NIR imaging become available for nearby 
galaxies. Our experiment with 2D masses that include NIR information thus confirms that the 
use of integrated photometry increases the random error in associated measurements for whole 
populations, while introducing a non-trivial bias as well. In the worst case, OD and 2D masses 
differ by ~0.25 dex when NIR data are involved, and ~0.20 dex otherwise. 


Fig. 12 can also be compared to Z09 and SS15, the only prior studies of the spatially-resolved 
approach for measuring A/*. Z09 found that, on a galaxy-by-galaxy basis, A/^^ could exceed 
by as much as ~40-65% (0.15-0.22 dex), depending on the MLCRs and dataset employed; 
the most egregious disagreements occured when they modelled optical-to-NIR SEDs. Similarly, 
SS15 found that, in the most extreme cases, measurements of missed ~25, 50, and 20% of the 
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stellar mass in E/SO, spiral, and Irr galaxy types. While we find offsets of similar magnitudes to 
Z09, they only oeeur in the sense of Ml^ > when we eonsider optieal information alone. With 
NIR fluxes ineluded, our 2D method typically yields smaller masses compared to OD measure¬ 
ments. A contributing factor might be that the BC03 models tend to predict optical-NIR MLCRs 
of shallower slopes than the CB07 models (Fig.[^, thereby leading to milder M*/Lgradients. This 
discrepancy raises further awareness of the need for tests of MLCRs and SPS models on resolved 
stellar populations; the biased distribution of points in Fig. [^pertaining to NIR-based predictions 
strongly hints at a disconnect between SPS models across the optical-NIR domain. However, it 
should be noted that SS15 measured (by as much as 25%) for nearly one-fifth of their 

full sample, with these galaxies being distinguished by their low-to-medium sSFRs. Recalling that 
our sample is drawn from a cluster environment and that cluster galaxies tend to exhibit lower 
sSFRs, relative to their field counterparts, the discrepancies we have uncovered may be genuine. 

Aside from some notable outliers, Z09 also argued in favour of a correlation between and 

galaxy type, whereby later-types have larger ratios. They surmised that late-type galaxies, on aver¬ 
age, contain more star formation and dust than early-types, hence the integrated light of the former 
should more often be biased to lower M^/L. A closer inspection of their Fig. 13 calls some as¬ 
pects of their claim into question, however. First, the authors do not present any statistical tests for 
the existence of a correlation (e.g. Spearman’s coefficient) and, second, the correlation disappears 
when NIR fluxes are included in their modelling. Our own results suggest no such correlation ex¬ 
ists for any combination of the SHIVir bands. Instead, our use of a more representative sub-sample 
of galaxies than Z09’s has demonstrated that the integrated light of all galaxy types, not just late- 
type spirals, can suffer from significant biases. SS15 also showed (their Fig. 3) that a galxy’s sSFR 
appears to be a better diagnotic of these biases. The fact that we recover a dearth of galaxies with 
may be a consequence of this relation and the fact that our sample consists entirely of 
cluster galaxies. A 2D analysis for a complete sample of galaxies is ultimately needed to establish 
the dependence of 


With a larger, more representative sample of galaxies, we confirm the core results of Z09 and 
SS 15 that MI^ can differ substantially for local galaxies. The degree of scatter seen in Fig. 12 
is comparable to that in Fig. 11 for all datasets considered, except gi SEDs; Ml^ scatter to a larger 
degree than in this case. These findings appear consistent with the suggestion that radial vari¬ 
ations in colour are the predominant cause of differences between spatially-resolved and integrated 
mass measurements, rather than azimuthal or stochastic (i.e. pixel-by-pixel) variations. However, a 
larger sample of is needed to further test this point. Until such time, our recommendation for 
the use of Ml^, or depend on the particular application. For accuracies <0.1 dex and/or 
small samples, we advocate a spatially-resolved approach (1D/2D), while statistical corrections for 
the use of integrated light may be invoked in all other cases (see also SS15). 




5. Conclusions 


Drawing on the literature, SPS models, and a representative sample of local galaxies, we have 
quantified the random + systematic uncertainties in stellar mass estimates derived from MLCRs. 
In particular, we have searched for possible biases in M* due to the use of MLCRs (as opposed to 
SED fitting), prior assumptions invoked in building MLCRs, scope of SED sampling, and mod¬ 
elling integrated versus spatially-resolved data. Both MLCRs and SED fitting are able to recover 
the masses of mock galaxies to within ~0.2 dex of scatter (even when limited to just g- and 
z-band photometry), and yield equivalent results when applied to real galaxies. We find that the 
priors dominate the systematics, contributing at most ~0.3 or 0.6 dex of uncertainty, depending on 
whether optical or optical-i-NIR photometry is used. Modelling integrated photometry is the next 
most important effect. In so doing, 0.06-0.07 dex of scatter is introduced in M* measurements, 
while an additional bias of -1-0.05-0.07 dex results from including NIR data. The latter effect may 
be an indication that SPS models lack consistency throughout the optical-NIR range. Prom the 
MLCR perspective, while the use of a single colour (e.g. g-i) to constrain M^/L leads to bi¬ 
ases on the order of ~10%, increased sampling of galaxy SEDs beyond a few filters spanning the 
optical-NIR range only serves to reduce the statistical errors in the results. 

Our investigation demonstrate that the MLCR and SED fitting methods yield comparable 
results, as expected since their assessment of M* both depend on SPS models. Therefore, improve¬ 
ments to the MLCR method should come from extensive analyses of resolved stellar populations in 
galaxies, whereby fundamental parameters are most accurately constrained via colour-magnitude 
diagrams. In the meantime, we advocate the use of the MLCR method for applications focussed on 
exceptionally large datasets and/or complex parameter spaces. Barring future evidence to the con¬ 
trary, we only recommend the use of a spatially-resolved approach to M* measurements for cases 
involving either small numbers of galaxies or if accuracies to better than 0.1 dex are required. 
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A. Appendix 

Table [^provides linear least-squares fits to the MLCRs ereated on the basis of the BC03 and 
FSPS SPS models (see Figs.|7]j^. Fits are provided for eaeh of the filter and eolour eombinations 
aeeessed by our dataset. The eoeffieients in this table follow the generie equation: \og{M^/L)x = 
mx X (colour) + bx- We omit the CB07 MLCRs as those ean be found in Table B1 of Z09. Table 3 
shows the Pearson eorrelation eoeffieient for eaeh of these fits. In all eases the eoeffieient is greater 
than 0.55, and ~70% of the time it is greater than 0.85. 
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Table 2. Coefficients of linear fits for the BC03 and FSPS MLCRs. 


Colour 

(1) 

nig 

(2) 

bg 

(3) 

nif 

(4) 

br 

(5) 

nii 

(6) 

bi 

(7) 

niz 

(8) 

bz 

(9) 


bn 

(11) 

BC03: 

g-r 

2.029 

-0.984 

1.629 

-0.792 

1.438 

-0.771 

1.306 

-0.796 

0.980 

-0.920 

g-i 

1.379 

-1.067 

1.110 

-0.861 

0.979 

-0.831 

0.886 

-0.848 

0.656 

-0.950 

g-z 

1.116 

-1.132 

0.900 

-0.916 

0.793 

-0.878 

0.716 

-0.888 

0.521 

-0.968 

g-H 

0.713 

-1.070 

0.577 

-0.870 

0.507 

-0.834 

0.454 

-0.842 

0.313 

-0.902 

r—i 

4.107 

-1.170 

3.325 

-0.952 

2.925 

-0.908 

2.634 

-0.912 

1.892 

-0.977 

r-z 

2.322 

-1.211 

1.883 

-0.987 

1.655 

-0.937 

1.483 

-0.935 

1.038 

-0.975 

r-H 

1.000 

-0.988 

0.814 

-0.809 

0.713 

-0.778 

0.634 

-0.786 

0.414 

-0.833 

i-z 

5.164 

-1.212 

4.201 

-0.991 

3.683 

-0.939 

3.283 

-0.931 

2.210 

-0.947 

i-H 

1.257 

-0.869 

1.024 

-0.713 

0.895 

-0.693 

0.792 

-0.706 

0.495 

-0.761 

z-H 

1.615 

-0.729 

1.316 

-0.600 

1.150 

-0.593 

1.015 

-0.616 

0.615 

-0.692 

FSPS: 

g-r 

1.897 

-0.811 

1.497 

-0.647 

1.281 

-0.602 

1.102 

-0.583 

0.672 

-0.605 

g-i 

1.231 

-0.805 

0.973 

-0.644 

0.831 

-0.597 

0.713 

-0.576 

0.426 

-0.592 

g-z 

0.942 

-0.764 

0.744 

-0.612 

0.634 

-0.568 

0.542 

-0.548 

0.316 

-0.565 

g-H 

0.591 

-0.655 

0.468 

-0.527 

0.398 

-0.494 

0.339 

-0.482 

0.191 

-0.515 

r—i 

3.374 

-0.745 

2.675 

-0.600 

2.275 

-0.556 

1.940 

-0.537 

1.120 

-0.554 

r-z 

1.795 

-0.670 

1.421 

-0.539 

1.206 

-0.502 

1.021 

-0.487 

0.570 

-0.513 

r-H 

0.824 

-0.542 

0.654 

-0.439 

0.555 

-0.418 

0.469 

-0.414 

0.254 

-0.463 

i-z 

3.709 

-0.550 

2.933 

-0.443 

2.484 

-0.419 

2.084 

-0.411 

1.112 

-0.457 

i-H 

1.073 

-0.460 

0.852 

-0.375 

0.722 

-0.362 

0.608 

-0.365 

0.322 

-0.430 

z-H 

1.493 

-0.414 

1.188 

-0.339 

1.008 

-0.333 

0.849 

-0.341 

0.449 

-0.417 
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Table 3. Pearson correlation coefficients for the BC03 and FSPS MLCRs. 


Colour 

(1) 

rg 

(2) 

rr 

(3) 

n 

(4) 

rz 

(5) 

rn 

(6) 

BC03: 

g-r 

0.98 

0.97 

0.96 

0.96 

0.92 

g-i 

0.99 

0.98 

0.97 

0.97 

0.91 

g-z 

0.98 

0.98 

0.97 

0.96 

0.89 

g-H 

0.95 

0.95 

0.94 

0.93 

0.82 

r—i 

0.98 

0.97 

0.97 

0.96 

0.88 

r-z 

0.96 

0.96 

0.95 

0.94 

0.84 

r-H 

0.90 

0.90 

0.89 

0.87 

0.72 

i-Z 

0.92 

0.92 

0.91 

0.90 

0.77 

i-H 

0.85 

0.85 

0.84 

0.82 

0.65 

z-H 

0.81 

0.82 

0.80 

0.78 

0.60 

FSPS: 

g-r 

0.95 

0.93 

0.91 

0.88 

0.74 

g-i 

0.96 

0.93 

0.91 

0.88 

0.73 

g-z 

0.95 

0.93 

0.91 

0.87 

0.70 

g-H 

0.94 

0.92 

0.90 

0.86 

0.66 

r—i 

0.94 

0.92 

0.90 

0.86 

0.69 

r-z 

0.93 

0.91 

0.88 

0.84 

0.65 

r-H 

0.91 

0.90 

0.87 

0.83 

0.62 

i-Z 

0.90 

0.88 

0.85 

0.80 

0.59 

i-H 

0.90 

0.88 

0.85 

0.81 

0.59 

z-H 

0.89 

0.88 

0.85 

0.81 

0.59 



